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Abstract. A multi- wavelength spatial and temporal analysis of solar high energy 
electrons is conducted using the August 20, 2002 flare of an unusually flat (71 — 1.8) 
hard X-ray spectrum. The flare is studied using RHESSI, Ha, radio, TRACE, and 
MDI observations with advanced methods and techniques never previously applied in 
the solar flare context. A new method to account for X-ray Compton backscattering 
in the photosphere (photospheric albedo) has been used to deduce the primary 
X-ray flare spectra. The mean electron flux distribution has been analysed using 
both forward fitting and model independent inversion methods of spectral analysis. 
We show that the contribution of the photospheric albedo to the photon spectrum 
modifies the calculated mean electron flux distribution, mainly at energies below 
~100 keV. The positions of the Ha emission and hard X-ray sources with respect 
to the current-free extrapolation of the MDI photospheric magnetic field and the 
characteristics of the radio emission provide evidence of the closed geometry of the 
magnetic field structure and the flare process in low altitude magnetic loops. In 
agreement with the predictions of some solar flare models, the hard X-ray sources 
are located on the external edges of the Ha emission and show chromospheric plasma 
heated by the non-thermal electrons. The fast changes of Ha intensities are located 
not only inside the hard X-ray sources, as expected if they are the signatures of the 
chromospheric response to the electron bombardment, but also away from them. 
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1. Introduction 



Analyses of the impulsive phase of solar flares detected in various elec- 
tromagnetic wavelength ranges are important for understanding the 
acceleration and propagation of particles. The accelerated electrons 
reveal their presence e.g. in hard X-ray and radio emissions, and cause 
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rapid heating of the chromosphere and subsequent emission in the 
hydrogen Ha line. 

One of the key questions is the form of the accelerated electron 
distributions. Their determination from the hard X-ray spectra depends 
not only on the emission mechanism but also on the processes affecting 
the propagation of the particles. As suggested by Brown, Emslie, and 
Kontar (2003), it is the mean electron flux distribution F{E) depend- 
ing only on the bremsstrahlung cross section that should be used as 
a reference for both observational hard X-ray spectral analyses and 
theoretical models of electron acceleration and propagation. The mean 
electron flux distribution can be determined either by inverting the 
photon spectra (Johns and Lin, 1992; Plana et al., 2003) or by a forward 
fitting method usually assuming a power-law form of F(E) (Holman 
et al., 2003). Observed hard X-ray spectra also contain a contribution 
of photons backscattered in the photosphere and this will modify the 
calculated F{E) (Bai and Ramaty, 1978; Alexander and Brown, 2002). 

Besides the X-ray emission, numerous types of radio emission in 
the metric and decimetric range are usually observed during flares. 
They provide further information about the acceleration and propa- 
gation of the electron beams and plasma parameters in the emission 
region (Karlicky, 1997). 

According to standard solar flare models (e.g. Dennis and Schwartz, 
1989), the accelerated electrons provide one of the mechanisms for the 
flare energy transport from the release site to the lower atmospheric 
layers. As the accelerated electrons propagate along the magnetic field 
lines toward the photosphere, they lose their kinetic energy mainly 
via Coulomb collisions (Brown, 1971) in the lower corona and chro- 
mospheric layers. The chromospheric response to the beam energy 
deposition determines the characteristics of optical and UV emission 
from the flare loops. Numerical models of chromospheric response to 
pulsed electron beam heating (e.g. Canfield and Gaylcy, 1987; Heinzel, 
1991) predict the time correlation of hard X-ray and Ha emission re- 
cently analysed e.g. by Trottet et al. (2000) and Wang et al. (2000). 
Their results show that hard X-rays and Ha intensities in some flare 
kernels exhibit time correlations in the time range from subseconds 
to ~10 s. 

Some flare models (e.g. Heyvaerts, Priest, and Rust, 1977; Cargill 

and Priest, 1983) and observations (Czaykowska et al., 1999) suggest 
that Ha kernels are located between the upflows of the beam heated 
plasma and the downflows of the cool plasma in the loops disconnected 
from the reconnection site. The leading edges of the Ha emission close 
to the newly reconnected loops arc supposed to be also heated by 
the accelerated particles, which are detectable e.g. as the hard X-ray 
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sources. Comparison of the spatial distribution of Yohkoh hard X-ray 
sources and Ha flare kernels was done by e.g. Asai et al. (2002). They 
found that many Ha kernels brighten successively during the evolution 
of the flare ribbons, but only a few radiation sources were seen in the 
hard X-ray images. They suggest that this discrepancy may be the 
result of the low dynamic range capability of the Yohkoh Hard X-ray 
Telescope. 

In this paper, we analyse the impulsive phase of the August 20, 2002 
solar flare, which is characterised by a very flat photon spectrum at the 
X-ray burst maximum (spectral index 71 = 1.8). Similar flat spectra 
have been reported, e.g. by Nitta, Dennis, and Kiplinger (1990) and 
Farnik, Hudson, and Watanabe (1997). The event presented here is 
the first such fiat spectrum spectrum observed with the high energy 
resolution of the Reuven Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI) (Lin et al., 2002). We show the influence of photons 
backscattered in the photosphere on the determination of the mean 
electron flux distribution and compare results obtained by inversion 
and forward fltting methods. We also study the spatial correlation of 
Ha emission and its time changes with hard X-ray sources. In section 
2 we describe the global behaviour of the event with regard to hard 
X-ray and radio fluxes. Sections 3 and 4 present the RHESSI, Ha and 
magnetic field observations, methods, and results. Finally, the results 
are discussed in section 5 and the conclusions are given in section 6. 



2. Global event description 

The analysed flare is well suited for a multi-wavelength study of the 
electron propagation because its observations provide a set of simulta- 
neous data of comparable time and spatial resolution in the wavelengths 
directly related to the accelerated particles: RHESSI hard X-ray spec- 
tra and images. Ha images (Kanzelhohc observatory), and radio and 
microwave fluxes (Bern, Ondfejov, and Ziirich facilities). 

The flare was detected by GOES satellites on August 20, 2002. It 
started at 08:25 UT and reached its maximum at 08:26 UT as an M3.4 
flare. The IB Ha flare was reported in NOAA AR 0069 at S10W38, 
starting at 08:25 UT, peaking at 08:26 UT, and ending at 08:37 UT^. 

A global overview of the flare evolution in dm/m radio waves and 
RHESSI X-rays is shown in Figure 1. The RHESSI flare was first de- 
tected at 08:24:32 UT with the increase in flux in the 12 - 25 keV energy 
band. Then, with some time delays (at ~ 08:24:45 UT), bursts in higher 

^ http : //www . sec . noaa . gov/weekly/pdf 2002/prf 1408 . pdf 
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08:24:30 08:25:00 08:25:30 08:26:00 08:26:30 

Stort Time (20-Aug-02 08:24:20) 



Figure 1, Time evolution of radio and hard X-ray fiuxos during the August 20, 
2002 flare from observations made with Ziirich and Ondfejov radiospectrometers 
and RHESSI. Narrowband emission (NbE), type III, and Illd bursts are annotated 
on the Ondfejov and Ziirich data plots. For display purposes the RHESSI count 
flux at the energy band 300 - 500 keV was multiplied by a factor of 2, radio fluxes 
were shifted with a step of log A = 0.2 from the fiuxos at 4299 MHz (Ondi^ejov) and 
988 MHz (Zurich). Attenuators states AO, Al, and A3 and the time interval in which 
the derivative of Ha intensity was analysed are indicated by the lines at the top of 
the RHESSI data plot. The abrupt changes in RHESSI flux in the 12 - 25 keV energy 
band at ~ 08:25:17 UT and ~ 08:25:41 UT are caused by the attenuator changes 
(AO — » Al and Al A3, respectively) at these times. 
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and higher energy bands followed. During the starting phase of the X- 
ray emission at 08:24:30 - 08:25:08 UT, no radio bursts were observed 
in the dm/m range. The first weak and narrowband emission (NbE) in 
the dm-range was registered at 08:25:08 UT in the 2000 - 2200 MHz 
frequency range, followed by the type III burst below 500 MHz at 
08:25:16 UT. 

The most energetic part of this flare occurred during the time inter- 
val 08:25:15 - 08:25:50 UT, when significant hard X-ray emission was 
detected up to 7 MeV (see also data from the SONG instrument on 
board CORONAS-F; Myagkova, priv. comm.) and microwave emission 
up to 89.4 GHz. The turn-over frequency between optically thick and 
thin part of the radio emission was shifted to values between 19.6 
and 50 GHz. At that time the radio emission in the dm-range also 
reached the maximum and a very fast drifting (~ -4000 MHz s~^) 
and short-lasting (~ 0.1 s) type Hid burst (relativistic type HI burst 
- see Poquerusse, 1994) was observed in the 400 - 1000 MHz range at 
08:25:31.8 UT. In the decimetric range (2000 - 4500 MHz) this phase 
was characterised by broadband pulsations. 



3. RHESSI data analysis 

RHESSI X-ray data were analysed with the RHESSI Data Analysis 
Software (Hurford et al., 2002; Schwartz et al., 2002). During the flare, 
different attenuators were automatically placed in front of all detectors 
to absorb soft X-rays and reduce pulse pile-up. The attenuator states 
AO, Al, and A3 (no attenuator, thin attenuator, and both thick and 
thin attenuator, respectively) are indicated on Figure 1. Corrections for 
the effects of the different attenuator states are available for all times 
except when the attenuators are moving, a period of ~ 1 s. 

3.1. Imaging 

Image reconstruction was performed with both the CLEAN and the 
PIXON algorithms using data from the front segments of the detectors. 
Detector 7 was used only for the images in the energy bands above 7 keV 
due to its energy threshold of about 7 keV. Detector 1 and 2 mainly 
added noise to the images indicating that there was no significant source 
structure with an angular scale finer than ~ 4 arcsecond. Therefore, 
these detectors were excluded. 

The CLEAN algorithm was used for obtaining morphology and time 
evolution of X-ray sources, see sections 4.1 and 4.3, Figures 5, 6, and 9. 
Since there was no source component as large as or larger than FWHM 
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of the detector 8 (107 arcsec), the imaging fields of view in all CLEAN 
images were set smaller than this FWHM and conseqTicntly detectors 
8 and 9 were also not used in the CLEAN image reconstruction. 

The PIXON algorithm (Puetter and Piha 1994; Metcalf et al., 1996) 
is known to suppress spurious sources and to have high photometric 
accuracy (Metcalf et al., 1996; Alexander and Metcalf, 1997). We have 
used the PIXON imaging technique for imaging spectroscopy at the 
burst maximum, see section 3.3. 

3.2. Spectroscopy 

We analysed data summed over the eight front segments. Detector 2 
was excluded since its energy threshold is about 20 keV and its energy 
resolution is about 10 keV. The data were corrected for pulse pilc-up 
and decimation. Full 2D detector response matrices, each corresponding 
to the applied attenuator state, were used to convert input photon 
fluxes to count rates. The response matrix accounts for the efficiency 
and resolution of the detectors, the absorption by the attenuators, 
grids, and all other material above the detectors, and all other known 
instrumental effects (Smith et al., 2002). 

The time bin for RHESSI spectra was 2 s (« a half rotation), which 
ensures that the rapid modulation produced by the grids does not have 
any distorting effect on the spectra. The RHESSI spectra were analysed 
in the time interval 08:25:10 - 08:26:00 UT and were fitted from 8 keV 
up to typically 400 keV, depending on the signal-to-noise ratio. 

RHESSI data of this flare were contaminated by an electron precip- 
itation^ starting at 08:16 UT and detected up to 300 keV. Therefore, a 
special care was taken to estimate the background due to the electron 
precipitation in the 12 - 300 keV energy range while standard tech- 
niques were used at other energies. In the 12-40 keV energy range, the 
background was estimated using count rates from the onboard charged 
particle detector whose lower threshold (50 keV) is generally triggered 
during electron precipitation events (Smith et al., 2002). Analysing 
several precipitation events which occured on other orbits two days 
before and after the flare, we have found a good correlation between 
the particle count rate and the front count rate in this energy range. 
Therefore, count rates from the particle detector can serve as a proxy 
for the bremsstrahlung of the precipitating electrons as detected by 
germanium detectors in this energy band. 

At higher energies the germanium detector count rates during the 
precipitation event are generally dominated by the bremsstrahlung that 

^ hsi_20020820_071620_rate.png and hsi_20020820_071620_part_rate .png 
at http : / /hessi . ssl . berkeley . edu/hessidata/metadata/ 2002/08/20/ 
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originates from the electrons interacting in the Earth's atmosphere and 
generally does not follow the time history of the particle count rate 
(Smith, priv. comm.). Having taken this into account, the precipitation 
component in the 40 - 300 keV energy range was estimated with a cubic 
fit to the shoulders of the precipitation event detected before and after 
the flare. 

The uncertainty of the estimated precipitation component in the 
12 - 300 keV energy range was estimated to be < 50%. Despite this 
high value, any residual precipitation component in the background 
subtracted count spectra in the 12 - 300 kcV energy range is < 10% at 
the edges of the analysed time interval and decreases to several percent 
with the increasing flare flux. Thus, the precipitation event does not 
affect our results in the analysed time and energy range. 

3.3. Flat photon spectrum at the flare maximum 

In order to evaluate the power-law indices of the photon spectra and to 
make comparisons with previous works, we adopted a standard photon 
spectrum model / consisting of a double power-law non-thermal com- 
ponent /lit (e.g. Lin and Schwartz, 1987) and a Maxwellian isothermal 
component Ith- I = ^th + -^nt- 



where e, eb, 7i, 72, and a are the photon energy, the break energy in 
keV, the power-law index below and above eb, and a is for normal- 
isation, respectively. We used the forward fitting method and fitted 
the background subtracted X-ray spectra with the model spectrum I. 
The model spectrum was convolved with the full 2D detector response 
matrix to compute the model count spectrum. Parameters of the best 
fit model were found by searching for a minimum Xmin between the 
measured and the model count rate spectrum. 

The lowest value of the power-law index 71 was found at the burst 
maximum as seen in > 25keV X-rays, at 08:25:22 UT, when it dropped 
to 1.8. Note that this value lies outside of the range of common values 
of photon power-law index, (2,7), as reported e.g. by Dennis (1985). 
The flatness of the spectrum does not result from the contamination by 
a generally flat spectrum of the precipitation event because the residual 
precipitation component in the background-subtracted count spectrum 
was estimated to be < 1% in the whole energy range at 08:25:22 UT. 
Hence, we believe that the flat spectrum is of solar origin. 

This conclusion is further supported by the results from the imaging 
spectroscopy using the PIXON imaging method (Hurford et al., 2002). 




(1) 
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RHESSI PIXON images are constructed by finding the simplest model 
which is consistent with the counts observed during the time inter- 
val. The model used included a source image and a component that 
accounts for the counts unmodulated by the grids. This additional 
component prevents the unmodulated counts from being included in 
the source image. Thus, the PIXON image is not contaminated by this 
background and the PIXON image reconstruction provides an indepen- 
dent method of determining the spatially integrated flare spectrum. We 
reconstructed the PIXON images in several energy bands for 4.141 s 
(~ 1 spin period) centred at the burst maximum (08:25:22 UT) and 
determined a PIXON photon spectrum from the total photon fluxes 
above the 10% contour level of the corresponding PIXON image. Since 
the uncertainty in the spectra from the image reconstruction is ~10% 
(Hmford, priv. comm.), and the PIXON photon spectrum agrees within 
uncertainties of the same order with the flat photon spectrum up to 
300 keV (i.e. the highest energy for which the sources were detectable), 
we may conclude that the background subtraction procedure, see sec- 
tion 3.2, resulted in a solar count spectrum not significantly distorted 
by the precipitation component. 

3.4. Correction for photospheric albedo and its 

IMPLICATIONS on MEAN ELECTRON FLUX DISTRIBUTION 

Primary photons with energies larger than 15 keV emitted downwards 
have a high probability of being refiected due to Compton backscat- 
tering in the photosphere. These albedo photons then modify both the 
spectrum and the intensity of the observed hard X-rays. According to 
Bai and Ramaty (1978), the reflectivity of a power-law photon spectrum 
from the isotropic source, which is assumed here, is a function of the 
photon energy and significantly depends on the photon spectral index 
7. The reflectivity peaks at ~ 35keV and increases as 7 decreases. The 
total photon spectrum, that is the primary one plus the component 
composed of the reflected photons, is then characterised by a lower 
7 than that of the primary spectrum. For a spectral index 7 ~ 1.8, 
the reflected photons at 35 keV can be as high as 80% of the primary 
flux. The amount of backscattered photons generally depends on the 
anisotropy of the primary hard X-ray source and the pitch-angle distri- 
bution of the beam electrons. But, the dependency is rather weak due 
to strong angular smoothing effect of the bremsstrahlung cross section. 
A detailed study of beam electron distributions is out of scope of this 
paper. However, the assumption of isotropy of the hard X-ray source 
provides a lower limit for the photospheric albedo correction. 
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For this flare we studied the influence of the photospheric albedo 
correction on the mean electron flux distribution F(E). We used two 
approaches. First, F{E) was approximated by a double power-law func- 
tion 

— ) ^E^^^ Elow < E < -Ebreak 

I "^-^break ^ -^break <E< i^^high 

where E, -Eiowi -E'break) -E'high) ^i-, <^2, and A are the electron energy, 
the low-energy cutoff, the break energy, the high-energy cutoff in kcV, 
the power-law index below and above -Bbreak; and normalisation, re- 
spectively. The mean electron flux distribution represents the density 
weighted electron flux distribution in the source, and its determination 
is equivalent to the determination of the electron flux distribution as- 
suming thin-target bremsstrahlung /thin(e) (Brown, Emslie and Kontar, 
2003) 

^high 

/thin(6) = ^ / nE)Q{^,E)<^E , (3) 

where n is the mean intensity in the source volume V , R = 1 AU is the 
distance of the emitting source, and Q{e, E) is the isotropic cross section 
for bremsstrahlung (Haug, 1997). Using the forward fltting method 
described in section 3.3, the background subtracted photon spectra 
were fltted with a Maxwellian isothermal component /th plus the thin- 
target bremsstrahlung /thin of the electron flux model F(E) and the 
best-flt values of thermal plasma and electron flux parameters (the flt 
provides nVF{E)) were obtained. .Bhigh was kept flxed at 7 MeV since 
setting the parameter free did not significantly improve the fits. 

Second, the mean electron flux distribution function was derived by a 
regularised inversion procedure (Plana et al., 2003; Kontar et al., 2004, 
2005a) using the photon spectra (determined with the forward fitting 
method) as input. Electron energies were determined in the energy 
range typically 10 - 400 keV, the upper energy depending on the signal- 
to-noise ratio. 



3.4.1. Photospheric albedo 

Background subtracted photon spectra were corrected for photospheric 
albedo using the method of an angle dependent Green's function for X- 
ray Compton backscattering (Kontar et al., 2005b), where the Green's 
function represents the probability density of backscattering of a pho- 
ton of initial energy cq into the observer's direction with energy e. This 
approach is independent of the spectral characteristics of the primary 
photon spectra and allows any form of the primary photon spectra to 
be deduced, not just a power-law spectrum as was done e.g. by Bai and 
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Figure 2. Influence of the photospheric albedo on the photon and mean electron 
flux spectrum at 08:25:22 UT (intergration time 2s). Top panel shows the total 

(crosses) and primary photon spectrum ( ) . Horizontal sizes of crosses correspond 

to widths of energy bands, vertical sizes represent ±lcr statistical plus systematic 
uncertainties. Best fit to the primary photon flux is composed of the isothermal 
-^th (• • •) and thin-target /thin (— • — ) component. Bottom panel compares forward 

fitted fiVF{E) derived from the total ( ) and primary ( ) photon spectrum. 

Low-energy cutoffs and their standard errors corresponding to the total and primary 
photon spectrum are El°l;^^ = 44 ± 6 keV and Ef^^ = 20 ± 20 keV, respectively. See 
section 3.4.2 for details on parameter error estimation. 



Ramaty (1978), and Alexander and Brown (2002). Green's function for 
a given heliocentric angle (40 degrees for this flare) can be applied to 
the primary photon spectrum but as a computational convenience we 
modify the detector response matrix to account for the photospheric 
albedo. In this way it can be used straightforwardly in the forward 
fitting method (Kontar et al.'^) because the result of the photon model 

^ http : / /www . astro . gla . ac. uk/user s/ eduard/rhessi/albedo/ 
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multiplied by this new modified detector matrix is exactly the same as 
the result from the original detector matrix multiplied by the product 
of the albedo and the photon spectrum. 

In addition to photospheric albedo, the RHESSI X-ray spectra gen- 
erally contain a component from photons scattered in the Earth's at- 
mosphere. The geometry of RHESSI detectors is such that most of the 
photons scattered in Earth atmosphere must first pass through the rear 
segments (see Figure 1 in Smith et al., 2002)., where the photons are 
effectively absorbed (the stopping depth of Ge detectors for a photon 
of 30 keV is ~ 1 mm). Thus, only the photons coming from a forward 
semisphere can get directly into the front segments. Since RHESSI 
is pointed towards the Sun, the front segments are shielded by the 
rear ones from the Earth scattered photons all the time when RHESSI 
can see the Sun apart from being at the terminator. The studied flare 
was observed when RHESSI was far from the terminator. Therefore, 
a correction for the Earth albedo photons for actual zenith angle was 
neglected in the analysis. 

The primary photon spectrum (derived under the assumption of 
its isotropy) together with the total photon spectrum at 08:25:22 UT 
(burst maximum) is shown in the top panel of Figure 2. The photon 
spectra differ mainly below 100 keV, the primary spectrum is slightly 
steeper, = 1-9) ^-iid exhibits less flattening below 50 keV than the 
total photon spectrum. Such an albedo correction significantly affects 
nVF{E) mainly at electron energies below 100 kcV. The bottom panel 
of Figure 2 shows that a forward fitted nVF{E) consistent with the 
total photon spectrum has a high value of E\°'^^ , which would represent 
a gap in energy between the thermal and electrons producing hard X- 
rays. However, removing the contribution of the photospheric albedo 
photons results in a forward fitted nVF{E) with Ef^^ close to the 
photon energy where the thermal part of the photon spectrum joins the 
non-thermal part. Consequently, no reliable dip in nVF{E) is found. 

3.4.2. Comparison of forward fitting and inversion methods 
The mean electron fiux distribution nVF{E) corresponding to the pri- 
mary photon spectrum at 08:25:22 UT derived from both the forward 
fitting method and the regularised inversion procedure are shown in 
the top panel of Figure 3. The ranges of parameter values were de- 
termined as a change — XminI ~ 8.2 while keeping the individual 
parameter fixed, where the level 8.2 encloses the region of joint con- 
fidence of 7 parameters with 68% confidence (Lampton, Margon and 
Bowyer, 1976; Press et al., 1992). Uncertainties of the forward fitted 
nVF{E) were determined by multiple generations of the forward fit 
parameters lying within the same level of joint confidence. zLl-a limits 
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of the regularised solution were found from multiple inversions of the 
primary photon flux data perturbed within corresponding statistical 
and systematic uncertainties. The inversion method was studied in 
detail in Kontar et al. (2004), where various modelled input spec- 
tra were analysed. Middle and bottom panels of Figure 3 represent 
normalised and cumulative residuals of the primary JP"" and modelled 
photon spectrum I = /th + -^thin- Normalised residuals are defined as 
r{e) = [IP''(e) — I(e)]/(T(e), where a{e) is the uncertainty in the primary 
photon spectrum and includes statistical plus systematic uncertainties. 
The error of the albedo correction method is ~ 6%, which is larger 
than (j(e) but not included in it. Cumulative residuals are defined as 
5(e) = l/NJ2fn, where is the iV*^ energy bin used in fitting, and 
can be used to assess clustering of residuals in certain energy ranges 
(Plana et al., 2003). Normalised residuals of the best fit are limited 
to a 3-c7(e) interval. Cumulative residuals show minor clustering in the 
photon energy range 20 - 80keV but they are well below ±l-a limits of 
random walk {a = \/\fN), see full lines in bottom panel of Figure 3. 
Both normalised and cumulative residuals also indicate that assumed 
systematic uncertainties (5% in each energy bin) are overestimated. 
Adopting a lower vahic of 3% results in residuals distributed as expected 
for Gaussian statistics with ~ 1/3 of them above the 1-a level, but it 
does not modify hVF(E) significantly. 

The best forward fit parameters for the primary spectrum are: emis- 
sion measure EM = (1 it Q4) x 10"^^ cm^^, temperature T = (1.3 it 
0.2) keV, ^low = (20 ± lO) keV, = 0.4 ± 0.2, 62 = 3.2 ± ^b.eak = 
(240 ± tg)keV, nVF = (17 ± g;^) x 10^^ electrons cm-2s-\ where F 
is F{E) integrated from -Eiow to -©high- Several parameters {EM, E'low; 
nVF) have significantly different, by a factor of 5, estimates of high 
and low a. This behaviour is mainly due to the low value of 5i because 
a decrease in Ei^-^ and consequent increase in nVF does not result in 
significant change in the modelled primary photon spectrum. The pri- 
mary photon spectrum is therefore, within uncertainties, consistent also 
with a mean electron flux distribution function without a low-energy 
cutoff as derived from the forward fitting method. Large uncertainties 
in EM are caused by the small number of data points in the thermal 
component of the spectrum (see top panel of Figure 2) . 

The mean electron flux distribution deduced from the forward fit- 
ting agrees within uncertainties with the regularised solution up to 
~ 200 keV. The discrepancy at higher energies is caused by different 
high-energy cutoffs in the forward fitting and inversion methods. The 
regularised solution was cut at 0.5 MeV whereas the double power-law 
F{E) extended to -Ehigh = 7 MeV. The large uncertainties between 20 
and 40 keV shown in Figure 3 result from the small number of electrons 
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Figure 3. Top panel shows comparison of nVF{E) deduced from forward fitting 

( ) and inversion (- - -) method for primary photon spectrum at 08:25:22 UT. 

Shaded area represents il-cr Umits on the regularised solution, thin full lines corre- 
spond to 68% confidence strip for the forward fit solution. Middle and bottom panels 
show normalised and cumulative residuals of the primary and modelled photon 
spectrum for the forward fitting method. Thin lines in the bottom panel represent 
±l-cr limits for a random walk process. 



in this energy range and the fact that the photons are primarily from 
higher energy electrons. A similar argument can be used to explain the 
deviation between the forward fit and regularised solutions above 100 
keV. The high value of i^high used in the forward-fit solution introduces 
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Figure 4- Time evolution of RHESSI X-ray count flux (1-s time resolution) in five 
energy channels, the power-law index 5i , and the low-energy cutoff i?iow of the mean 
electron flux distribution nVF{E). Attenuator states AO, Al, and A3 are indicated 
by lines at the top of RHESSI data plot. For display purposes the RHESSI count 
flux at the 300 - 500 keV band was scaled by a factor of 2. 

a contribution to the photon flux above 100 keV from the higher energy 
electrons not included in the regularised electron spectrum. 

In the case of the total photon spectrum, i.e. without the albedo 
correction, the regularised solution indicates a dip within 1-a limits in 
the energies below ~ 30 keV where the forward fitted solution resulted 
in a gap above the thermal component in nVF(E), see Figure 2. 

3.4.3. Time evolution 

Figure 4 shows the time evolution of derived parameters 6i and -Ejo^ 
of the mean electron flux distribution corresponding to the primary 
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photon spectra. The power-law index Si decreases to its minimum of 
0.4 ± 0.2 at the burst maximum, 08:25:22 UT, and its time history 
generally shows a soft-hard-soft pattern between the photon flux and 
hardness of the spectra (e.g. Fletcher and Hudson, 2002). Low values 
of Si indicate a flat electron flux spectrum and this is confirmed by the 
regularised inversion method. The low-energy cutoff -Blow lies within or 
close to the energies where the thermal component dominates nVF(E) 
distribution. Standard errors are of the order of the E\ow values. There- 
fore, no reliable dips or gaps in nVF{E) throughout the time interval 
could be established. Our analysis shows that the low-energy cutoffs 
above the energies of the thermal component are likely to be the result 
of the contribution of photospheric albedo photons and do not exist in 
real mean electron spectra. 



4. Ha data analysis and magnetic field topology 

To study the influence of the accelerated electrons on the Ha emission, 
we used full-disk Ha images covering the whole interval of the flare. The 
images were acquired using a filter with 0.07 nm FWHM at Kanzelhohe 
Solar Observatory in Austria. The time resolution of the data is 3 s and 
the spatial resolution is 2.2 arcsec/pixel. The position of the solar centre 
and the spatial resolution were determined using the routines provided 
by A. Veronig. Details on the method can be found in Veronig et al. 
(2000). 

The intensities from the flare site were scaled by the averaged quiet- 
Sun area intensity Iq for every image separately. The brightest parts 
of the flare area beginning at 08:25:20 UT are distorted by the satu- 
ration of detected intensities. The saturated pixels are displayed with 
horizontal lines in the Ha images (Figures 5, 6, and 9). 

The data set was further spatially coaligned using a technique devel- 
oped for TRACE and EIT data (Gallagher^). The cross-correlation of 
the Ha data set was done in the area of the nearby sunspot with respect 
to the chosen Ha image. That sunspot area was not affected by the flare. 
Therefore, we assumed no time change of the Ha intensities during the 
analysed time interval in that area. The positional uncertainty of the 
Ha images was estimated to be of the order of 1 pixel (2.2 arcsec). 

4.1. Spatial evolution 

The evolution of the hard X-ray sources overlayed with the Ha emission 
and MDI magnetic neutral lines is displayed in Figures 5 and 6. The 

^ http : / /hesperia . gsf c . nasa . gov/~ptg/trace-align/ index . html 
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Figure 5. Time and spatial evolution of Ha intensities (grey scale) and X-ray emis- 
sion (full line contours) in the 7-12 keV (first column) and 12 - 25 keV (second 
column) energy ranges. Contours correspond to 50, 70, and 90 % of maximum of each 
RHESSI CLEAN imago. There were no detectable X-ray sources in the 12 - 25 keV 
energy range at 08:24:39 UT (see the image in the top right corner). MDI magnetic 
neutral lines are plotted as dashed lines. The saturated Ha pixels are displayed with 
horizontal lines. 
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25 - 40 keV 



40 - 70 keV 
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X (arcsec) 

Figure 6. Time and spatial evolution of Ha intensities (grey scale) and X-ray emis- 
sion (full line contours) in the 25 - 40 keV (first column) and 40 - 70 kcV (second 
column) energy ranges. Contours correspond to 50, 70, and 90 % of ma^ximum of each 
RHESSI CLEAN image. There were no detectable X-ray sources in the 40 - 70 keV 
energy range at 08:24:48 UT (see the image in the top right corner). MDI magnetic 
neutral lines are plotted as dashed lines. The saturated Ha pixels are displayed with 
horizontal lines. 



RHESSI flux was accumulated for 1 spin period (4.141 s). The Ha 

images correspond to the middle of the RHESSI imaging time range. 
Due to the saturation of the Ha intensities, the spatial evolution in the 
flare area is shown to 08:25:20 UT only. 

At the very beginning of the flare (at 08:24:39 UT, Figure 5), a single 
X-ray source is detected. Then it splits into two (at 08:24:48 UT, Fig- 
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ures 5 and 6) and appears to be moving eastwards. Later at 08:25:14 UT, 
the X-ray emission starts to shift southwards and the soTirccs spht even 
into three components detectable in the energy range 40 - 70 keV. There 
is a difference in the positions of hard X-ray sources at lower and higher 
energy ranges detected at the same time. The X-ray sources at higher 
energies tend to be located more to the south than the lower energy 
X-ray sources (compare the positions of the sources in different energy 
bands at 08:25:14 UT and 08:25:20 UT in Figures 5 and 6). Most of X- 
ray source maxima at energies above 12 keV are not spatially coincident 
with the brightest parts of the Ha emission. Rather they tend to be 
located near the north and eastern edges of the bright Ha kernels. The 
Ha kernels form a semi-circle whose inner boundary is located near the 
magnetic neutral line. 

4.2. Magnetic topology 

To understand better the magnetic field topology near the flare site, 
the extrapolation of magnetic field in the current-free approximation 
(Demoulin et al., 1997) was made using the MDI magnetogram ob- 
served at 08:12:30 UT (see Figure 7). 




560 575 590 arcsec 



Figure 7. Magnetic field lines, extrapolated in the current-free approximation at 
08:12:30 UT (full thick lines), revealing the quadrupolar topology near the flare site. 
The isolines of the north and south magnetic polarities (± 10, ± 100, and ± 400 G) 
are expressed by thin full and dashed lines, respectively. A cut through the region 
of the neutral magnetic line indicated by the arrow is shown in the sketch in the top 
right corner (here sizes of loops are not to scale) . 
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The most interesting structure, indicated by an arrow in Figure 7, 
is located along the magnetic neutral line, which is delineated by the 
above mentioned semi-circular Ha ribbon. Namely, this structure has 
a quadrupolar character in the sense of papers by Uchida (1996) and 
Hirose and Uchida (2002), i.e. low altitude loops in the close vicinity 
of the neutral line are between magnetic field lines with the opposite 
orientation and all is covered by extended overlying loops. 

The preflare activity in the 1550 A TRACE band was observed 
directly in the centre of the quadrupolar structure (i.e. above the mag- 
netic neutral line) at 08:12:20 UT. However, no 3 - 6 kcV X-ray emission 
was detected by RHESSI from this flare site at this time. 

Comparing Figures 5, 6 and 7, it can be seen that the flare started 
at the north end of this quadrupolar structure. Then the flare kernels 
appeared at the north-cast boundary of the strTicturc and the flare ends 
at the southern extremity of the quadrupolar structure. We also note 
that the Ha ribbon and X-ray sources are shifted ~ 4 arcsec from the 
quadrupolar configuration of the extrapolated magnetic field lines. 

4.3. Time derivative of Ha intensities 

We assumed that the locations of the optical flare that are directly 
affected by the electron beams might be recognised as areas of abrupt 
changes of the Ha intensity, /h«, on a time scale of the order of seconds. 
This assumption is supported e.g. by Trottct et al. (2000), who found 
correlation between hard X-ray flux and Ha intensities on time scales 
from seconds to tens of seconds. Therefore, the time derivative of lua 
may provide a technique to find such areas. 

In order to suppress the effects of seeing and/or incorrect positions. 
Ha intensities were averaged over a 3 x 3 pixel area and the resulting 
lightcurves were smoothed by a 3-point box car (9 s time resolution). 
Such smoothing was necessary to distinguish between the real changes 
and those which occurred randomly throughout the whole flare area. 
The effects of the averaging and smoothing are shown in Figure 8. 
Due to the saturation of Ha data beginning at 08:25:20 UT, the time 
derivative of Ha intensity, dlna/dt, was studied in the time interval 
08:24:33 - 08:25:20 UT (see also the label in Figure 1). The maximum 
of d/Ha/dt was 0.097 IqS-^ A threshold, dlua/dt = 0.03IqS-^ was 
estimated by fitting d/na/dt distribution outside of the flare area with 
the Gauss distribution. The threshold corresponds to 10 standard de- 
viations from the mean of the fit and was adopted as the level of the 
reliable value of d/na/dt not affected by noise or seeing. 

The locations of d/na/dt > 0.03 IqS"-*- are displayed on Figure 9 
together with the RHESSI CLEAN images in the 25 - 40 keV and 40 - 
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Figure 8. Averaging and smoothing efTects of /ho at the location [560,-280] arcsec 
reveahng the largest derivative of Ihc (see Figure 9). Thin lines show observed time 
histories at 3x3 pixel area, () correspond to the average over that area, and the thick 
line is the resulting time history smoothed by 3-point box car. Intensities are scaled 
by the quiet Sun intensity 7q. 



70 keV energy bands. X-ray sourees in the 40 - 70 keV energy band 
are not shown in Figure 9(a) since no source at this energy band was 
detectable at that time. 

These energy bands were chosen for comparison with d/no/dt for 
two reasons. First, they were above the energies where the thermal 
component dominates X-ray emission, i.e. > 12keV. Secondly, the 
electrons bombarding the chromosphere need to have an energy greater 
than 25 kcV to penetrate into the layers of Ha formation, e.g. in the 
Fl flare model by Machado et al. (1980). The RHESSI CLEAN images 
were reconstructed from data accumulated in 1 spin period. The centre 
of each RHESSI time interval corresponds to the observational times 
of Ha and d/Ho/dt. 

The sites of the fast Ha intensity changes located in the northern 
part of the flare lie on the boundaries or inside of X-ray sources in the 
25 - 40 keV and/or 40 - 70 keV energy range (see Figures 9(b) and 
(c)). However, we have found other areas that reveal abrupt changes 
in Ha intensity but are located well away from X-ray emission (e.g. 
Figure 9(a)). Some of these areas even precede the position of X-ray 
emission as it moves to the south in time (Figures 9(c) and (d)). The 
difference in the locations of d/na/dt and the X-ray sources cannot be 
explained by a spatial correlation with the X-ray sources at lower or 
higher energies. There is no significant shift of the X-ray sources in the 
energy ranges from 3-7 keV to 7 - 12 keV or from 40 - 70 keV to 
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Figure 9. Locations of time derivative of Ha intensities (full line contours) overlayed 
on Ha images (grey scale). Contours correspond to 0.04 IqS~^ and 0.09 IqS~^ of 
d/Ha/dt. RHESSI CLEAN images at 25 - 40 keV and 40 - 70 keV are plotted 
as dotted and dashed contours at 50 and 90 % of maximum of each image. The 
saturated Ha pixels are displayed with horizontal lines. 

higher energy bands; the sources at energies from 7 to 25 keV are not 
located closer to the d/Ha/dt (see Figures 5 and 6) either. 

The fact that at 08:25:20 UT there are no locations of d/na/dt > 
0.03 Iq s~^ along the whole flare area may be the result of the saturation 
of the Ha data at this time (areas covered by horizontal lines in Fig- 
ure 9). Although the described analysis technique does not eliminate 
the saturated pixels, they only affect the adjacent pixels and cause a 
decrease in d/na/dt. 



5. Discussion 

The photon spectrum at the burst maximum of the August 20, 2002 
flare is, to our knowledge, the flattest (71 = 1.8) yet detected by 
RHESSI with its high energy resolution. To eliminate any distortion of 
the spectra, the background subtraction also included the estimation 
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of the electron precipitation component, which has been done for the 
RHESSI data for the first time. The solar origin of such an unusu- 
ally flat spectrum was confirmed by a comparison with the spectrum 
constructed from the RHESSI PIXON images. 

Using a new method of the photospheric albedo correction indepen- 
dent of the spectral form of the primary spectrum, we have proved 
the importance of the photospheric albedo correction in the case of the 
analysed flare. We have shown that in this case of fiat X-ray spectra, the 
photospheric albedo significantly affects the determination of the mean 
electron flux distribution nVF{E), mainly the value of the low-energy 
cutoff. The contamination of the flare spectra by the photospheric 
albedo photons may result in dips in fiVF{E) indicated both by the for- 
ward fitting and inversion methods. Such dips could be misinterpreted 
as spectral features inherent to the acceleration/propagation processes 
of the energetic particles in the solar flares. 

Flares of flat hard X-ray spectra (7 < 2) were previously analysed 
e.g. by Nitta, Dennis, and Kiplinger (1990) and Farnik, Hudson, and 
Watanabe (1997). All of these flares have weak thermal components 
(GOES class M3 or lower) but intense non-thermal components up to 
several hundred keV with a flat part extending down to 10 keV. As 
was shown here, such photon spectra may correspond to a flat electron 
flux distribution with the low-energy cutoff close to the energy of the 
thermal electrons. 

For the flares with spectral index around 3 it is the low-energy 
cutoff which determines the total energy content in the flare. This 
happens when the largest contribution comes from the low energy 
particles. However, for the flares with flat spectra < 1 it is not 
the case. Assuming the collisional thick-target scenario (Brown 1971), 
one can estimate the energy flux of the non-thermal electrons injected 
into the plasma. In this model, the injected electron flux distribu- 
tion !F{Eq) [electrons s~^ kev"-*^] is related to F{E) of Equation (2) as 
(Brown and Emslie, 1988) 

J'iE^) ~ F{Eo)E^^ (d In F{E)/d lnE + 1)^^^^ ~ ^(^o)^o~' • (4) 

Relative change of the non-thermal energy flux ^Ei^^ [keV s~^ ] can be 
expressed 

^E{JdE,^^^{Ei^/E,,^)-'^ 5biow=/ HEo)EodEo. (5) 

Adopting the spectral index 5i = 0.4 and the low-energy cutoffs, E[^^ = 
10 keV, -Eiow = 30 keV, consistent with the inferred flat spectrum, the 
ratio is ~ 2. Thus, due to the flatness of nVF{E), the low values of the 
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low-energy cutoff do not cause a significant increase in the non-thermal 
energy flux. 

The flat distribution function of non-thermal electrons, especially at 
low energies, implies that its derivative in the energy as well as in the 
velocity space is close to zero or even positive. Such a distribution 
function can be formed directly in the acceleration region (e.g. by 
the acceleration in the direct electric field - Holman, 1985) and/or 
modified by propagation effect, particle collisions in dense plasma, and 
by interactions of accelerated electrons with plasma waves (Melnik 
and Kontar, 2003). If the plasma waves are generated, then the radio 
emission can be observed at high frequencies, as in our case. However, 
detailed conditions in the radio source needed to confirm this idea are 
not known to us. 

Our analysis of the characteristics of the hard X-ray sources and the 
Ha and radio emission revealed several facts about the flare process. 
First, the presented flare was weak in the metric range during the whole 
event, except for the short-lasting rclativistic type HI burst (radio 
manifestation of rclativistic electrons; Klassen, Karlicky, and Mann, 
2003). It resembles the March 6, 1989 7-ray flare, which started as 
silent in radio waves below 800 MHz. Similarly to Rieger, Treumann, 
and Karlicky (1999), we explain the lack of metric radio emission by 
closed geometry of the magnetic structure and a low altitude of the 
flare process. Additionally, the high value of the turn-over frequency, 
19.6 - 50 GHz, of the radio emission indicates a strong magnetic held 
or a large depth of the emitting source (Benka and Holman, 1992), 
thus supporting this conclusion. Furthermore, the coalignment of the 
Ha and hard X-ray images with the magnetic field extrapolation shows 
that both the Ha kernels and the hard X-ray sources are located at the 
footpoints of the low altitude magnetic loops which are the part of the 
quadrupolar structure (Figure 7). 

Second, our results concerning the time and spatial correlation of 
the hard X-ray sources and Ha emission are partly in agreement with 
the prediction of solar fiare models. We have found that most of the 
X-ray emission maxima, mainly those at the energies above 12 keV 
(corresponding to the X-ray emission of the non-thermal electrons), 
tend to be located on the northern and eastern edges of the Ha emission 
(see Figures 5 and 6). This is consistent with the expectation that X- 
ray sources should be located on the external edges of Ha ribbons and 
similar to those reported by Czaykowska et al. (1999). 

Flare models also suggest that fast changes of Ha intensities might 
be located at positions of X-ray sources, i.e. at the areas which are 
bombarded by the suprathermal electrons. In some cases, this sugges- 
tion is confirmed by our analysis. However, at some times the positions 
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of the fast changes of Ha intensities do not spatially coincide with the 
X-ray sources (see southern locations of dlna/dt in the bottom row 
of Figure 9). Wc propose that such Ha brightcnings might be caused 
either by thermal conduction from a hot plasma located above, or by 
soft X-ray irradiation. Using the lowest energy ranges of RHESSI (3 - 
7 keV), wc tried to find such a hot plasma at these places, but without 
any success. However, we must point out that the time resolution of 
diHa/di was 9 s. Therefore to confirm the existence of these areas, more 
analyses and better quality data are needed. 

Finally, the time evolution of the positions of the hard X-ray sources 
may indicate the movement of the reconnection site. This can be seen 
e.g. at 08:25:14 UT (Figures 5 and 6) when the hard X-ray sources of 
different energies are distributed along the Ha ribbon in such a way 
that the sources of higher energy arc located further to the south- 
east direction. Such difference in the positions could be expected in a 
single flare loop when it would correspond to the height distribution 
of the hard X-ray sources along the loop as found in Aschwanden, 
Brown, and Kontar (2002). However, this is not the case here because 
the extrapolation of the magnetic field shows no magnetic loop along 
the Ha ribbon. The displacement of the Ha ribbon and X-ray sources 
from the quadrupolar configuration is not due to the positional errors 
in the Ha, MDI, and RHESSI images, which factor of 2 and 

4 smaller, respectively. The offset could be caused by neglecting the 
electric currents in the current-free extrapolation. Although the linear 
force-free extrapolations show similar results, we cannot exclude the 
presence of localised currents which deform the magnetic field structure. 



6. Conclusions 



We have demonstrated that the influence of the photospheric albedo 
on solar flare X-ray spectra must be considered in order to correctly 
assess the models of electron acceleration and propagation and gen- 
eration of hard X-ray emission. The inferred flat mean electron flux 
distribution could be the result of the acceleration process and/or mod- 
ified by particle-wave interactions. This explanation, however, remains 
speculative. The accumulated information from the multi-wavelength 
observations of this flare shows that the flare process took place in the 
low altitude magnetic loops. Our analysis of Ha emission may indicate 
that the response of the lower atmosphere to the flare energy release is 
not restricted to the sites of propagation of the accelerated electrons. 
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